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We construct twisted instanton solutions of CP" models. Generically a charge-A: instanton splits into k{n + 1) well-separated and 
almost static constituents carrying fractional topological charges and being ordered along the noncompact direction. The locations, 
sizes and charges of the constituents are related to the moduli parameters of the instantons. We sketch how solutions with fractional 
total charge can be obtained. We also calculate the fermionic zero modes with quasi-periodic boundary conditions in the background 
of twisted instantons for minimally and supersymmetrically coupled fermions. The zero modes are tracers for the constituents and 
show a characteristic hopping. The analytical findings are compared to results extracted from Monte-Carlo generated and cooled 
configurations of the corresponding lattice models. Analytical and numerical results are in full agreement and it is demonstrated 
that the fermionic zero modes are excellent filters for constituents hidden in fluctuating lattice configurations. 

•Key words: cooling, instanton, overlap operator, supersymmetry, topology, zero mode 
'PACS: ll.27.H-d, ll.lO.Wx 
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1. Introduction 

■Nonlinear sigma models in two dimensions have long been used 
3S testing ground for strongly coupled gauge theories They 
are scale invariant on the classical level and asymptotically free 
3t the quantum level. The ubiquitous CP" models possess reg- 
ular instanton solutions, the topological charges of which yield 
lower BPS-bounds on the action, they have a chiral anomaly 
when coupled to fermions, generate a dynamical mass by non- 
perturbative effects at zero temperature and a thermal mass oc 
g^T at finite temperature. They have numerous interesting ap- 
plications to condensed matter physics (for a review see 12(]) and 
also have been used to study the sphaleron induced fermion- 
number violation at high temperature fst]. 

In the present work we consider CP" models at finite tem- 
perature, i.e., Euclidean models with imaginary time having pe- 
riod P = l/fefiT. These models possess instanton solutions 
,with finite action and the dimension of the moduli space in a 
given instanton sector depends on the topology of the Euclidean 
space-time. For example, on the two-torus the charge-fc instan- 
tons of CP" depend on as many collective parameters as the in- 
stantons of CP"+^ with one charge less 101 ■ In the present work 
we do not compactify space such that space-time is a cylinder. 

For suitable field variables the selfduality equation for CP" 
instantons reduces to Cauchy-Riemann conditions such that all 
instantons are known explicitly for the plane, cylinder and to- 
rus. On the plane they are given by rational functions of the 
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complex coordinate z and on the cylinder by suitable periodic 
generalizations thereof jSee below. 

In a previous work ||5|] one of us introduced the twisted 0(3) 
model (which is equivalent to the CP^ model) and showed that 
generically the unit charged instantons in this model dissoci- 
ate into two fractional charged constituents, sometimes called 
'instanton quarks'. Again there is a close analogy to the cor- 
responding situation in Yang-Mills theories, where instantons 
with nontrivial holonomy along the compact direction of a four- 
dimensional cylinder possess magnetic monopoles as consti- 
tuents SHiii. 

We extend the work in jst] in several directions. First we 
construct the k{n 1) constituents of CP" instantons with 
charge k and twisted boundary conditions and relate their po- 
sitions, sizes and fractional charges to the collective parame- 
ters of the instantons. Then we calculate and analyze the zero 
modes of the Dirac operator for minimally coupled fermions 
with quasi-periodic boundary conditions in the background of 
the twisted instantons. We show that the zero modes can be 
used as tracers for the instanton constituents: they are local- 
ized to the latter, to which constituent depends on the boundary 
condition. Again, this has close analogies in four-dimensional 
Yang-Mills theories with 1 (or 2) compact dimensions 



CP" spaces admit a Kahler structure such that the two- 
dimensional CP" models admit a supersymmetric extension 
with two supersymmetries. These models contain 4-fermi in- 
teractions and the Dirac operator is given by the linearized field 
equation for the fermions. We calculate the zero modes of this 
operator. There exists always one zero mode with squared am- 
plitude being proportional to the action density of the twisted 
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instanton. 

We supplement our analytic studies by numerical simula- 
tions. With known algorithms we produce typical field config- 
urations for various CP" models with twisted boundary condi- 
tions. Then we apply standard lattice cooling techniques to ex- 
tract the instantons and their constituents from a given (thermal- 
ized) configuration. Again we find that an instanton of charge 
k consists of k{n + 1) constituents. How the corresponding in- 
stanton constituents in Yang-Mills theories emerge in the pro- 
cess of cooling/smearing has been studied in ifisl [l4l [Tsi [3, 
TtI 181. Next we compute and analyze the zero modes of the 
overlap Dirac operator We find good agreement between our 
numerical and analytical results. 

Our results demonstrate that CP" models and Yang-Mills 
theories share one more common feature: in both models the 
twisted instantons generically split into well-separated and al- 
most static constituents and in both models the fermionic zero 
modes trace these constituents. In SU (N) Yang-Mills theories 
the number of constituents is given by k and the rank whereas 
for CP" models it is given by k and by n. On the lattice the 
fermionic zero modes are excellent filters for the constituents. 
Even without much cooling the zero modes detect the con- 
stituents of the fully cooled configurations. 

2. The CP" model in the continuum and on the lattice 



The two-dimensional CP" model | 2^ can be formulated in 
terms of a complex (n + l)-vector u = {uq, . . . , subject 
to the constraint u'^u = 1. The Euclidean action is given by 



S 



It is invariant under local U{1) gauge transformations 

uj[x) ^ e'^(^) Uj[x), A^,{x) ^ A^{x) + a^A(x), 
as well as global transformations 

Uj{x) t-^ Ujiui{x) 



(1) 



(2) 



(3) 



with a constant matrix W G [/(n + 1). The gauge field can 
be eliminated from the action by using its algebraic equation of 
motion. 



A., 



-iv)d^u. 



(4) 



The integer- valued instanton number, 

Q = Jd^xq{x) with (7(2;) = ^e^^a^A^(a;), (5) 

can be interpreted as the quantized magnetic flux in a fictitious 
third dimension. At infinity u must approach a pure gauge. 



u{x) 



c and Q is just the winding number of the map 
at infinity, an element of the first homotopy group 
of [/(I). 

Configurations minimizing the action in 5 > 4:TT\Q\/g'^ are 
called instantons. They fulfill first order self-duality equations. 



The most general instanton solution with instanton number Q 
fc G N can be written in homogeneous coordinates Vj as 



u-i{x) 



\v{z)\-' 



i = 0, . . . , n, 



(6) 



with {vj } a set of polynomials of the complex coordinate z = 
xi + 1x2 with no common root and maximum degree k. The 
topological charge density of an instanton configuration then 
reads 

q{x) ^ ^A\n\viz)\\ (7) 

Lattice formulation 

For the bosonic model the lattice regularization can be obtained 
as described in 11211 12211 . After introducing the matrix-valued 
gauge invariant field 



P{x) = u{x)u''{x), 



(8) 



which projects onto the one-dimensional subspace spanned by 
u, one finds 

tr [d^P d^P] = 2d^u^d^u + 2 {u^d^^uf = g^C. (9) 

This equation, valid for the model defined on a continuous space- 
time, is discretized naively with the forward derivative, S^P 1-^ 



P 



Px, such that 



tr [d^P d^P] ^ 2d - 2 ^ tr [P,P,+^] . 



(10) 



Therefore, the action, up to an irrelevant additive constant, takes 
the form 



2 ^ 2 ^ 



(11) 



The simulations of the lattice models have been performed with 
the help of an overrelaxation algorithm ||22i] . In addition, to in- 
vestigate the topological properties, we cooled the lattice con- 
figurations |23[] . For a given configuration one cooling step con- 
sists of minimizing the action locally on a randomly chosen site 
X. This is achieved by constructing = J2fj,iPx+fi + Px-fi) 
and replacing by the eigenvector corresponding to the largest 
eigenvalue of Q^. A cooling sweep corresponds to one cooling 
step per lattice site on average. Using this procedure the instan- 
ton constituents naturally emerge from the locally fluctuating 
fields. 

For the topological charge on the lattice we used the geo- 
metric definition in |0] leading to an integer-valued instanton 
number This definition and the chosen lattice action are suf- 
ficient for the analysis of global topological properties in the 
vicinity of classical configurations. Thus, we are not affected 
by the improper scaling behavior of the dynamical CP" mod- 
els withn < 2 111- 
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Figure 1: (Color online.) Logarithm of the topological density for the 1-instanton solution of the CP^ model (see and )19t ) with symmetric constituents, 
/^i = ^2 — Ml = 1 ~ ^^2 = 1/3 (cut off below e~^). The parameters Ai are chosen such that the constituents are localized according to i22\ from left to right at 
(ai, 02, as) = (-5, 0, 5), (-5, 1, 4), (-5, 7, -2) (first line) and (-1, 0.5, 0.5), (0, 0, 0), (3, -1, -2) (second line). Note that the xi-range has been changed 
in the lower right panel. 



3. Instantons at finite temperature 

For a quantum system at inverse temperature f3 we identify z ^ 
z + i/3. Since (3 is the only length scale in the problem we 
measure all lengths in units of /3. In particular the coordinates 
become dimensionless and we identify z ~ z + i. Periodic 



fc-instanton solutions ('calorons') are given by 11261127 1 



^per 



(z) = 6(0) + e2-~' + . . . + 6^=) . 



,27Tkz 



(12) 



By a global U{n + 1) symmetry transformation one can rotate 
Wpei- such that the constant (and per assumption non-vanishing) 
vector fe('^) G C"+^ points in the 0-direction, bj'") = b^^SjQ. 

The twisted model is only quasi-periodic in the imaginary 
time direction. This means that the components vj of v are 
periodic up to phases e-^'^'^J with /ij G [0, 1), i.e., the vectors 
V and u are periodic up to a diagonal element of the global 
symmetry U{n + 1). The U{n + l)-invariants like \v\ and 
and hence also q stay periodic. Without loss of generaUty we 
assume the phases to be ordered according to /xq < /ii < . . . < 

For the general solutions of the twisted model we consider 
the Fourier ansatz 



Hence the total topological charge in the twisted model can 
have a fractional part, whose values are restricted by the bound- 
ary conditions. By a global transformation we enforce Kmin to 
be taken on in the 0th component and by a (non-periodic) local 
transformation we further set yUo = and = 0, such that 
Q = Ktnax- According to Eq. these powers also govern the 
asymptotic values of the fundamental fields uj . 

In the following we will mainly analyze twisted instantons 
with integer-valued instanton number Q = k E N. They are 
obtained by K^ax = k taken on in the 0th component, i.e., the 



(fc) 



highest coefficient b^'''> points in the 0-direction, b 
Thus one can obtain the components Vj by multiplying each 
component Wper,j from ( fT2] i with cxp{2TTfijz), which yields 



v{z) = 51 Wper(z), 51 = diag (e 



(16) 



For n = 1 the known twisted unit charged instanton solu- 
tion isl can be recovered in terms of the gauge invariant field 



vo(z) 



^0 

(1) 



(17) 



We made use of /ig = and b\ = and denoted /ii by w. 



E 



(13) 



is) 

and demand the coefficients 6^ to be non-vanishing only for 
a finite range of s (for each component j). This is because the 
corresponding maximum and minimum of the powers 



max (s + /^j ) , 



mm 



3.1. One-instanton sector 

In order to explore the topological density of the instantons we 
first consider solutions with unit charge Q = k = 1. We mul- 
tiply w by a constant such that b'^^ = 1 and afterwards shift the 
Euclidean time X2 such that b'^^ becomes real and non-negative. 
For this choices the density |wp only depends on the absolute 
(s + /ij) (14) values Xj = \bj''^ \ with j = 0, 1, . . . , ti. If, in addition, we de- 



then yield a finite topological charge Q. According to (|7]l one 
has to compute the following surface integrals 



d.T2 5i In 



a;i — i-oo 
a;i— > — cx 



— ^max '^min 

G No + {mj - I ^ = 0, . 



(15) 



i}. 



fine An 



fep^) I and fin+i ~ 1, then it can be written in the 



condensed form 



\viz)\' 



71+1 

J2 \} e^-^'^i 



2A„+ie2''^i cos(27r.T2) 



(18) 



The corresponding topological charge density splits into n + 1 
constituents at most. For CP^ this is illustrated in Fig.[T|which 
shows h\q{x) for various choices of the parameters A^. 



For the general CP" models the occurrence of the con- 
stituents can be understood geometrically. To see this more 
clearly we write 



n+l 



2 6^^(^1) 008(277X2), 



i=0 



with 



Pi{xi)~ AirniXi + 21nAj;, 
p{xi) = 271X1 + In A„+i. 

In particular 

Pa{xi)= 0, 
p„+i(xi)= Attxi + 21nA„+i = 2p{xi] 



(19) 



(20) 



(21) 



We compare the graphs of these n + 3 linear functions, see 
Figs. |2}Slfor three examples in the CP^ model amounting to 
five exponential terms. 

The dominant contribution to |wp in ( fT9] l at a fixed point 
xi comes from the exponential term whose graph is above the 
fines defined by the other exponential terms. Hence In is 
piecewise linear in the direction xi up to exponentially small 
corrections that are maximal in transition regions, where the 
highest lying graphs intersect. 

Note that for a strictly linear In | u | ^ the topological density 
q oc A In would vanish exactly, whereas at cusps generated 
by intersections of linear parts the topological density would 
be a Dirac delta distribution (in 3+1 dimensional Yang-Mills 
theory a similar singular localization can be obtained in the far- 
field limit l28l ll ill ). As this is a good approximation to the 
actual In we conclude that the topological density of the 
twisted instantons splits into constituents localized at the inter- 
section points of the lines. Because of the ordering of the /i's, 
the slopes of the linear functions pi are ordered as well. Note 
that for xi < — l/(27r/ii) In Ai the term exp(po) dominates 
such that In |w| « on the left of all constituents. Correspond- 
ingly, \n\v\ « Attxi on the right of all constituents. 

We obtain the maximum number of constituents, if all con- 
secutive graphs intersect separately and above the rest of the 
graphs, respectively. More precisely said, the twisted instanton 
of CP" splits into + 1 constituents, if, and only if, ai <^ 
02 ^ • ■ • <C a„+i{3 whereas ai is the intersection point of the 
fines pi^i and pi, 



ln(Ai/Ai_i) 



2n (fii 



fJ-i 



l)' 



1, 



(22) 



i.e., in particular the twist parameters /i; are distinct, = /iq < 
^1 < • • • < /i„ < fin+i = 1- These positions are arbitrary 
provided the corresponding A-parameters are chosen according 
to 



In A, 



-2nY,{f^i 

1=1 



l-H-i)ai with /io = 0. 



(23) 



' Thereby we do not want to take the condition a; _ i <^ a; too literally. It is 
sufficient, if „ i is not close to , whereas the required distance is determined 
by the slopes ofp^ - pi_i andpi_i - Pi-2- 
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Figure 2: In and exponents pi and p as a function of xi, see Eqs. 11 9t - 
)20t . in the CP^ model for the case of (ai , 02, as) = (—5, 1, 4), which leads 
to three well-separated constituents (equivalent to 2nd example in Fig.^. 
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Figure 3: In | tip and exponents p; and p as a function of xi, for the case of 
(ai , a2 , as) = (—5, 7,-2), where the second and third constituent merged 
(equivalent to 3rd example in Fig.^. 
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Figure 4: In and exponents pi and p as a function of xi, for the case 
of (ai,a2,a3) = (3,-1,-2), where the time-dependent p-term becomes 
relevant (equivalent to 6th example in Fig.[T). 
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Let us discuss the case of the well-separated constituents. In 
the neighborhood of the intersection point of the lines pi-i 
and Pi we can approximate 

|2 



+Xt e 



2 4TTfliXi 



(24) 



for xi not too far from the constituent i, 

- (ai_i + aj) <xi <- (ai + a.^+i) , i = 1, 



,n+ 1. 



(25) 

For i = 1 the lower bound for xi is — oo and for i = n + 1 the 
upper bound is +cxd. The contribution to the topological density 
of the ith constituent is 



(26) 



cosh^ [2TT{fii - Mi-i) [xi - fli)] 

This shape is the same for all CP" models, cf. jst] for the CP^ 
case. The constituent decays exponentially with characteristic 
length 2TT{fii — fJ-i-i) (measured in units of (3) away from its 
position Qi. It has a fractional topological charge Qconst.i = 

— fii-i and these charges add up to 1 as they should. In 
terms of the linear graphs the fractional topological charge is 
proportional to the difference of slopes of the lines that meet 
(which also would give the amplitude of the delta distribution 
mentioned above), and the total charge is the sum of all slope 
differences, which indeed bend the graph from pq with slope 
to Pn+i with slope Att eventually. 

Neighboring constituents can merge adding up the fractional 
topological charges. This can be understood as 'pulling down' 
the line that connects the two constituents in the graph of In 1 1; [ ^ , 
in other words by choosing the corresponding parameter Ai 
small (cf. Figs. |3}Hll. 

Under which circumstances does the time-dependence of 
jwp contained in the last term of Eq. [19] (which is proportional 
to expp(xi)) play a role? Since the three graphs of pq, Pn+i 
and p intersect at the point (5,0) with 27ra = — lnA„+i we 
have 

p{xi) < max{po(2;i),Pri+i(a;i)} , (27) 

such that the time-dependent p-term can contribute to the sum 
in (T% only in the neighborhood of a. Furthermore, all other 
hnes have to lie below (a, 0) (if one of the other lines pi lies 
well above that point, then the topological density becomes to 
a good approximation static). As the slopes of the pi{xi), i ~ 
1, . . . , 77 are between and 1, these graphs are never dominant 
once they are below (a, 0). This means that only one transition 
point occurs. Hence time-dependence of the instanton appears 
iff the topological charge is concentrated in one lump (which 
can be thought of as all constituents merged, cf. Fig.|4]i. 

The case of non-distinct /^i's can be understood by consider- 
ing the limit fij f^j-i for some j's. Then the corresponding 
constituent becomes flatter and broader, in the limit it will be in- 
visible and 'massless' (i.e., without topological charge/action; 
this has been 'eaten' by the constituent j + 1). In this spirit 
we also recover the periodic solutiorH with /^o = Mi = • • = 



/i„ = 0: The resulting topological density then consists of only 
one constituent with unit charge, which can, but does not have 
to be time-dependent (depending on where the invisible, mass- 
less constituents are localized, see also the first row in Fig. 1 of 
Istl). This can be demonstrated by means of Fig. |2l i.e., based 
on the case of well-separated constituents: The limit is taken by 
sending the slopes of the graphs of pi and p2 to 0. If all posi- 
tions a-i are kept constant (by adjusting the A-parameters), then 
Pi, 2 (as functions) and the resulting topological density 
of the periodic solution is equivalent to the case of all con- 
stituents merged in the twisted model, cf. Fig. |4] and Fig. [T| 
(lower right panel). If the hmit is taken with all A-parameters 
kept constant (i.e., by sending the positions ai.2 to — oo), then 
Pi, 2 2 In Ai,2 can lie well above in the limit and the topo- 
logical density remains static, though it consists of only one 
unit charged constituent. 

Finally, we want to mention the possibility of generating 
solutions with topological charge less than 1 from these instan- 
tons. Technically one has to avoid the asymptotics ln|7;| w 
Attxi for large xi (we want to stay in the gauge where Kmin = 
and hence keep the asymptotics for smaU xi). Hence, if the cor- 



responding parameter A„ 



(1) 



is vanishing (in the Fourier 



ansatz there is no integer phase c^'^^), then the total topologi- 
cal charge of the configuration is less than 1. A phase with 
Kmax < 1 then gives the total topological charge (i.e., governs 
the asymptotics for large a;i). Also these configurations consist 
of constituents with the same formulae for locations, sizes and 
charges. The number of constituents varies from n down to 1, 
depending on how many of the remaining parameters A are zero 
(in the graph the corresponding lines and intersection points are 
missing). 

Interestingly, all these configurations have in common that 
the constituents in them are ordered along the noncompact di- 
rection. This has already been observed in jstl and substantiated 
by topological considerations. Here it can best be understood 
from Eq. (l22l i. The fractional charge of the ith constituent, 
/ii — 1 , is fixed by the twist in the boundary condition. These 
charges can be realized in isolation only if the ordering of their 
positions, oi <C ■ • • "C a„+i, applies. If some do not obey 
this ordering, then constituents emerge with the sum of the in- 
dividual fractions as their topological charge. In other words, 
'pulling a constituent through a neighboring one' results not in 
a different ordering but in joining the constituents to a bigger 
one, cf. Fig.[T](upper and lower right panels). 

Notice that by giving up our choice K,mm = we can rear- 
range the constituents cyclically; this can become relevant on 
the lattice, where xi is of course periodic. 

This ordering is of course related to the selfduality which 
dictates that all solutions are functions of xi + 1x2; antiself- 
dual solutions will have the opposite ordering. We therefore 
believe that this phenomenon is particular to 1h-1 Euclidean di- 
mensions. 



-In Yang-Mills theories this amounts to the Harrington-Shepard caloion 
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Figure 5: (Color online.) Logarithm of the topological density for the charge-2 instanton of CP^, with non-symmetric constituents, /Ji = 0.55, ^2—^*1 = 
0.15,1-/^2 = 0.3. The positions of the constituents from left to right are (ai, a2, 03, 04, as, ae) = (—10, —6, —2, 2, 6, 10), (—10, —4, —4, 2, 6, 10), 
(-6, -6, -6, 2, 6, 10) (first line) and (-10, -4, -2, 0, 6, 10), (-10, -2, -2, -2, 6, 10), (0, 0, 0, 0, 0, 0) (second line). 



3.2. k-instanton sector 

The general twisted instanton solution (fTZt . (fTSI l with integer- 
valued topological charge k corresponds to the norm 



\v{z)Y 



fc(n+l) 

E 

i=0 



k (fc-s)(n+l) 



'(^1) 



i=0 



X COS {2tISX2 + LPi+s{n+l) - Vi) , (28) 

where we introduced 

Pi{xi) = 47r/i,xi + 21nAi, 
Pi'^H^i) = 27r(2^,; + s)xi + In (AiAj+s(„+i)) (29) 

= ^ +P,;+s(ri-Hi)(a;i)] • 

We encoded the two indices of into one, i = s{n + 1) + j, 
and defined 



A,; 



(30) 



where s = 0, . . . , fc, j = 0, . . . , n. To arrive at (l28l l we trans- 
formed the constant vector 5'*^' in the 0-direction. Similarly as 
for the one-instanton solution we conclude that the constituents 
are localized at the transition points of the piecewise linear 
function [^(z)!^. The topological density thus splits into at most 
fc(n+ 1) constituents. Well-separated constituents are static and 
exponentially localized at a.i, i ~ 1, ■ • • , k{n + 1), given in an 
analogous manner as in the 1 -instanton case, cf. Eq. (l22l) . The 
constituents carry the fractional charge fii — and from the 
periodicity of the /i's in Eq. (l30l l follows that they are again 
ordered, see also Fig.|5] upper left panel. 
Again we have 



p[^\xi) < max{pj(xi),pj+s(n+i)(a;i)} , 



(31) 



(s) 

since the three graphs of pi, Pi+s{n+i) pl all intersect at 



a. = 



In {\+s{n+l) / ^i) 



2t:s 



(32) 



Therefore the time-dependent term containing pj'^' only con- 
tributes to the sum in ( l28b if the s{n+l) constituents at the posi- 
tions Qi+i, . . . , ai^s{n+i) merge to one constituent with integer 
charge s. The integer s thus determines the maximal frequency 
of the emerging constituent measured in units of the smallest 
possible frequency. We illustrated this behavior in Fig. |5] for 
2-instanton solutions of the CP^ model. Note that the freedom 

( s) 

ipi of choosing the complex phase of the parameters bj enters 
only as shifts in the time dependence, Eq. ( |28] ). 

3.3. Cooling of lattice data 

Charge-one instantons of the CP^ model containing up to the 
maximal number of three constituents are reproduced with a 
cooling of our lattice data. The simulations were performed on 
a 6 X 100 (temporal X spatial) lattice at coupling g^^ = 2. In 
the spatial direction periodic boundary conditions are imposed 
whereas in the temporal direction the vj are twisted with pre- 
scribed Hj . Then a particular configuration is cooled. During 
the cooling procedure configurations with = 1 and two or 
more well separated constituents are fairly stable even with this 
type of unimproved cooling (at least up to 10'' cooling sweeps). 
We also observe the typical annihilation of selfdual and anti- 
selfdual constituents. Only a small fraction of configurations 
is cooled to a state with three clearly separated constituents. 
One of these is shown in Fig.|6]at three different cooling stages. 
More often we end up with only two constituents. These results 
indicate that in a dynamical simulation topological objects with 
fractional charge (given by the twist parameters /ij) may be as 
relevant as they are in Yang-Mills theories. 

4. Zero modes of the Dirac operator 

4.1. Minimal coupling to fermions 

We extend the bosonic CP" model by introducing a massless 
Dirac fermion if}, for the time being minimally coupled, such 
that the action has the form 



S = 



A^x [{D^,u)^D^,u + iin'^D^^] 



(33) 




Figure 6: (Color online.) The cooling procedure is applied to a CP'^ configuration of the Monte-Carlo ensemble with twists fii = 0.15, ^2 — Ml =0.5, I — fJ.2 = 
0.35. Three stable constituents emerge after several cooling sweeps (10, 25, 500 from left to right). Here only the positive part of q is shown. 



We shall use the chiral representation of the 7-matrices for which 
7* = — i7^7^ is diagonal, 
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V 

1 0; ' 



7 



i^ 

-i or 



7* 



-1 
1 



(34) 



Splitting the Dirac spinor into chiral components ip ~ {ip, x)^ 
the Dirac equation in the background of a self-dual configura- 
tion splits into the Weyl equations 



(9 — u^duj if ~ \v\ d (^\v\ ^ if^ = 0, 
(i9 — v)du) X = B {\v\ x) ~ 0, 



(35) 



where d and d denote the derivatives with respect to the com- 
plex coordinates z = xi+ 1x2 and z = xi— \X2- It follows that 
the zero modes have the form 



ip{x) = f{z) \v\ and xi^) 



9{z) 



(36) 



with (anti-)holomorphic functions f{z) and g{z). Similarly as 
for the bosonic fields we impose quasi-periodic boundary con- 
ditions for the fermi field, 

il}{xi,X2 + l)^e^'''H{xi,X2) with Ce[0,l). (37) 

The functions /, g can be expanded in Fourier series. 



(38) 



s — — 00 



These modes are only square integrable on the cylinder if the 
coefficients a'^^\l3'''^^ and the twist parameter fulfill certain 
constraints. Recall that the asymptotic behavior of the general 
solution with charge Q = K^ax (we set Kmin = 0) is 



lim = 1 and lim |?;| (x e' 

Xi — ' — 00 Xi — !-00 



2tiQxi 



(39) 



Therefore there are no normalizable left-handed zero modes <p0 
The quantum number s of the right-handed zero modes is con- 
strained hy Q < s + C, < Q. 



^Vot anti-instantons with negative topological charge the left-handed modes 
become normalizable. 



ln|t;|^ 






47r(s + 
















1 .^ + c ■■■■■ 

















ai 



a2 



Xi 



Figure 7: Piecewise linear function In |t>| and linear function 47r(s + C)^l 
for the example of s + ^ G (/^i, ^2)- For x\ < a\ the slope of In is 
47r/Jo = 0, for ai < xi < 02 it is 47r/ii, for a2 < x\ it is 47r^2. The 
vertical distance between the graphs is minimal at 02 where the zero mode is 
localized. 



This immediately leads to the iyidex theorem for right-han- 
ded zero modes. For integer topological charge we have to dis- 
tinguish two cases for the fermionic phases C,: 



C-0 
Ce(0,l) 



Q ~ \ zero modes, 
Q zero modes. 



(40) 



For fractional topological charge we introduce the floor func- 
tion [ ■ J : M Z and the fractional part { ■ } : M -> [0, 1) such 
that Q = [(5J + {(3} and obtain 



C = o 

CG(0,{Q}) 

Ce[{Q},i) 



[QJ zero modes, 
YQ\ + 1 zero modes, 
YQ\ zero modes. 



(41) 



Let us further investigate the localization properties of the 
(right-handed) zero modes in the background of the instanton 
with integer charge Q = k. They have the explicit form 



^4Tv{s+Qxi-\n\v\' 



s = 0, . . . , fc - 1, (42) 



with \v\ from ( I28I I. It is helpful to first consider well-separated 

|2 



constituents for which In \ v 
piecewise linear. 

In \v{z)\^ « {Piixi) = ^TTHtXi + 21nAj I ai < xi < a^+i} , 

(43) 



becomes time-independent and 
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Figure 8: (Color online.) Minimally coupled fermionic zero modes (yellow/light) in the background of 1-instanton constituents (red/dark) of the CP'^ model with 
symmetric constituents, fj,i = fJ.2 — fJ-i =1 — ^2 = 1/3. In the first line the constituents are well-separated at (ai , 02, as) = (—5, 0, 5), in the second line from 
left to right they are at (—2.5, —2.5, 5), (—2.5, —2.5, 5), (0, 0, 0), i.e., two resp. three constituents have merged. The fermionic twist is ^ = 1/6, 1/4, 1/3 (first 
line) and 2/3, 5/6, 1/2 (second line). 



as described in Sec. 13.21 Clearly, the zero mode has maximal 
amplitude at points where In — 47r(s + C,)xi is minimal. At 
these xi the vertical distance between the graphs of the approx- 
imately piecewise linear function In and the linear function 
47r(s + C)xi is minimal. For s + C in the interval {fii-i,^i) 
the minimum is at a; 1 = Ui where the graphs of and pi 
intersect. The situation is depicted in Fig. |7] 

Hence, for generic values of ( the zero mode is localized 
at one constituent. The profile of the zero mode is symmetric 
about the constituent i for the particular value s + C, = + 



cosh [27r(/ii - ^i^i){xi - ai)]' 



(44) 



Interestingly, the profile is almost constant between the ith 
and {i + l)th constituent for ^ = /i^. These 'bridges' can be un- 
derstood by the fact that in this region the graphs of In and 
47r(s + C,)xi are parallel (up to exponentially small corrections) 
at these values of s and the fermionic phase C,. 

Altogether the zero modes walk along the ordered set of 
constituents when changing the fermionic phase With phases 
at the bounds C resp. C, = l(orC = {Q} for configura- 
tions with fractional charge) the zero modes become constants 
asymptotically. In other words, these zero modes have 'bridges' 
coming from ~oo resp. reaching out to +00 and hence are not 
normalizable. 

Similar arguments apply if two or more constituents merge. 
A few examples are given in Fig.|8] 

4.2. Zero modes on the lattice 

Using the overlap operator with quasi-periodic boundary con- 
ditions for the U{1) gauge field we are able to analyze its zero 
modes llsoll in the background given by the cooled lattice con- 
figurations. Even for moderately cooled configurations do the 
zero modes reflect the position of the fully cooled instanton 
constituents for specific boundary conditions (see Fig.|9]l. There- 
fore the lattice results are in full agreement with the analytical 



results and in addition the fermionic zero modes are excellent 
tracers for instanton constituents for mildly cooled configura- 
tions. 

4.3. Supersymmetric coupling to fermions 



The supersymmetric CP" model Oil 13211 contains n + 1 Dirac 
fermion fields ipj, j ~ 0, . . . , n, in addition to the complex 
scalar fields w, . Its action is 



S 



+ 



1 



where u and tp are constrained by 

u^u — 1, ■u'lp — ipu — 0. 



(45) 



(46) 



Introducing Weyl spinors, tp = {ip, x)^, the model is invariant 
under the on-shell JV = (2, 2) supersymmetry transformations 



Su = Slip - 82X1 

5ip ~ +2ieiDu - ei{(p<p)u + e2(x^)u, 
5x = -2\e2Du + £2(xx)w - £i(<^x)m, 
with the covariant derivative 

D = d — v)du and D ~ d — v)du 



(47) 



(48) 
£1 and 



and anticommuting parameters ei ^ satisfying £3 
e\ = £2- Both spinors ip and x have n + 1 components. 

The linearized Dirac equation in an external u-field splits 
into two Weyl equations. 



1 - uuf) Dp = and (1 - uu^) Dx = 0. 



(49) 



For an instanton background with u ~ v{z)/\v\ the Weyl equa- 
tions simplify to 



{1 - P„) d (\v\-' p) ^(1-P,)d{\v\x) = 0, 



(50) 




Figure 9: (Color online.) Fermionic zero modes of tlie overlap operator for the CP'^ model with twists = 0.15, /i2 - 
of 25 cooling sweeps (compare to Fig.|6] middle panel). The fermionic twist is = 0, 0.075, 0.15 (first line) and = 



-fj,i =0.5, l — fJ.2 = 0.35 after application 
0.4, 0.65, 0.825 (second line). 



where projects onto the holomorphic v{z), 

P.^uu^^^. (51) 

It follows that a left-handed solution reads 

ip{x) = \v\ fiz), (52) 

where /(z) is an arbitrary vector of anti-holomorphic functions 
orthogonal to v. None of these solutions is normalizable. With 
the help of PydPy = dP^ one shows that a right-handed solu- 
tion has the form 

x{x) = ^Al-Pv)9{z), (53) 

where g{z) is a vector of holomorphic functions. In order not 
to break supersymmetry g must fulfill the same boundary con- 
ditions as the instanton solution v. Therefore, the choice of 
fermionic twists is very limited here. Each function g can be 
constructed by linear combination of the basis elements {.g'^'*' } 
defined by 

g(^'^)(z)=e2-(«+''^)^e„ j = 0,...,n, s G Z, (54) 

where Cj is the unit vector pointing in direction j. For the cor- 
responding zero modes the squared norm is 

Ix^^'f = , ,4 El^H - (55) 

Normalizability of the zero mode in the fc-instanton background 
requu-es 

/0,l,...,/c forj = 0, 
s = < (56) 
[0, for j = 1, . . . ,n. 

In the case of well-separated instanton constituents we can write 

k(n+l) 
I'^^l ; mod {n+l)^i 




+ Mo 



fli a2 03 xi — > 

Figure 10; Logarithm of denominator and numerator of |x'^' |^ in i57\ . The 
zero mode has two maxima of equal amplitude at ai and 02. 

where we introduced x'*' = x'-''^-' for i = s{n + 1) + j. The 
linear functions pi{xi) are given in (|29] l. Again the maximum 
of Ix'^*-' I is easily found by considering the graphs of the linear 
functions 2pi{xi) and pi{xi) + An^iXi. In a logarithmic plot 
both the numerator and denominator of jx''' | are piecewise lin- 
ear For xi < ai the slope of the numerator is larger and for 
xi > Gi+i the slope of the denominator is larger. This is il- 
lustrated in Fig. [To] Simple geometric arguments about these 
graphs reveal, that the zero modes x^*-* with < z < k{n + 1) 
split into two constituents located at a; and Oi+i, which have 
the same amplitude, but decay with different lengths. The zero 
modes x'"' and x'*^'""*"^-" have only one maximum at ai and 
o-k{n+i)^ respectively. Some examples are plotted in Fig. (TT] 
The general right-handed zero mode has the form 

fe(n+l) 

Its (squared) norm splits into k{n+l) or less constituents. They 
have the same analytic form oc cosh~^ [27rQconst,i(a^i — ai)] 
and are located at the same positions ai as the instanton con- 
stituents. 




Figure 11: (Color online.) Supersymmetric coupled fermionic zero modes (yellowAight) in the background of 1-instanton constituents (red/dark) in the CP'^ model 
with twist parameters = 0.55, ^2 — /^l = 0.15, 1 — ^J,2 = 0.3. In the first line the zero modes x'"'. x'^^ ^nd x'^' (from left to right) with instanton locations 
(ai , a2, as) = (—5, 0, 5). In the second line the half-BPS state Xinsi ™th instanton locations (—5, 0, 5), (—5, 2, 3), (0, 0, 0) (from left to right). 



There exists always a particular zero mode, whose (squared) 
norm is proportional to the topological density 

1 I ,2 1 \v\^\dvf -Iv'fdvl^ 

9=:;- Ainu ^ (59) 

47r TT l^p 

Namely, since the squared norm of x in ( l53T l is 



,2 \v\'\g\' 



we obtain the exact relation 



\xix)\ =TTq{x) 



(60) 



(61) 



for the zero mode with g = dv, which means that /3i oc /i; A; in 
Eq. dSlll. 

The occurrence of this particular zero mode can also be un- 
derstood as follows: Any instanton background breaks half of 
the supersymmetry, namely the one generated by the parameter 
£2. If we transform the configuration tiinst = v{z)/ \v\, ipinst ~ 
with the broken symmetry then (JV'inst is inevitably a zero mode 
of the Dirac operator for the action is invariant under the super- 
symmetry transformation. This way we obtain a non-vanishing 
right-handed zero mode 

SXinst = ~2\£2Du OC jij (1 - P„) dv, (62) 

\v\ 

that is a zero mode with g oc dv. Except for irrelevant prefac- 
tors, this is exactly the zero mode whose squared norm is equal 
to the topological density of the instanton. 

5. Conclusions 

In the present work we constructed and analyzed the integer- 
charged instantons for twisted CP" models on a cylinder The 
twisted instantons with charge k support k{n + 1) constituents. 
If these constituents are well-separated then they become static 
lumps. The fractional charges and the shapes of the constituents 



topological profile are governed by the phases in the boundary 
condition (and the scale jj). The constituent positions are re- 
lated to the collective parameters of the twisted instanton and 
hence free up to the demand that for all constituents to be pre- 
sent their positions must be ordered. 

Neighboring constituents can merge adding up their char- 
ges. If at least n + 1 constituents merge then the resulting 
lump becomes time-dependent. For a composite object con- 
taining multiples of 71 + 1 constituents time-dependent terms 
with higher frequencies contribute, respectively. 

Our analytic findings are in complete agreement with the 
corresponding numerical ones. The latter were obtained by 
cooling lattice configurations of the twisted model with a non- 
vanishing topological charge. 

We determined all fermionic zero modes in the background 
of the twisted instantons. This has been achieved for minimally 
coupled fermions satisfying quasi-periodic boundary conditions 
in the Euclidean time direction. We found that, similarly as for 
gauge theories, the zero modes are localized at the positions of 
the constituents and that they may jump from one constituent to 
the neighboring one if the boundary conditions for the fermions 
are changed. Again we compared our analytical findings to nu- 
merical results. To that aim we determined the zero modes of 
the overlap Dirac operator for lattice configurations with differ- 
ent degrees of cooling. Again we find full agreement between 
our analytical and numerical results, in close analogy with the 
corresponding situation for SU (N) Yang-Mills theories. 

In the supersymmetric CP" model the Dirac fermions trans- 
form according to the fundamental representation of the global 
U{n + 1) symmetry group. The linearized field equations for 
the 71 + 1 fermion-flavors define a supersymmetric Dirac op- 
erator. We studied the square integrable zero modes of this 
operator and showed that they generically split into k{n + 1) 
constituents with maxima at the locations of the instanton con- 
stituents. There exists always a particular zero mode whose 
norm squared is equal to the topological charge density of the 
supporting instanton. This zero mode is generated by the half- 
broken supersymmetry. We did not elaborate on the contribu- 
tion of the constituents and zero modes to the central charge of 
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the (2,2) SUSY algebra. 

Our results are in close parallel to the corresponding find- 
ings in SU (N) gauge theories. But since twisted instantons, 
their constituents and the fermionic zero modes in CP" models 
are much simpler as in gauge theories our results may be use- 
ful to shed further light on the relevant degrees of freedom in 
strongly coupled models at finite temperature. The next natu- 
ral step would be to include quantum fluctuations about twisted 
instantons to study the quantum corrections to the constituent 
picture. 

In the SU (N) gauge theory there is a beautiful construc- 
tion of the constituents based on the Nahm transform fl. We 
believe that a similar construction, with Nahm transform as in- 
troduced in ||4i], could further simplify the construction of in- 
stanton constituents for twisted CP" models. 

Similar aspects of twisted CP" models (such as loop groups) 
are discussed in a simultaneous paper Issil . 
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